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We expose an analysis of the magnetic field distribution in the Abrikosov lattice of high-w; 
superconductors with d-wave pairing in the case where the critical field is mainly determined by the 
Pauli limit and the superfluid currents mainly come from the paramagnetic interaction of electron 
spins with the local magnetic field. The result found in frame of the generalized Clem variational 
approach is compatible with the recent observation that the form factor in CeCoIns increases with 
increasing field and then decreases at the approach of Hc2- (A. D. Bianchi et al., Science 319, 177 
(2008)). 
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Recent neutron scattering experiments performed on 
the heavy-fermion superconductor CeCoIns have re- 
vealed an unexpected behaviour of the vortex lattice 
(VL) form factoi'^"^ defined as the Fourier transform of 
the local magnetic field in the vortex lattice. The VL 
form factor of the type II superconductors is usually a de- 
creasing function of the magnetic field.'^ On the contrary, 
the VL form factor in CeCoIns was found to increase 
with increasing mag net ic field and then to fall down at 
the approach of i?c2 CeCoIns is a tetragonal, d-wave 
pairing superconductor with a large GL parameter, and 
with the highest critical temperature (T^ = 2.3 K) among 
all the heavy fermion compounds.^^lt^ It has already gen- 
erated great interest caused by the signs of the existence 
of the FFLO phase for a magnetic field parallel to the ab 
plane (and possibly to the c-axis),!^ and by the presence 
of an interval where th e sup erconducting/normal phase 
transition is first order.'^Efl' An explanation of the form 
factor behaviour has been proposed,^- which was based 
on a numerical processing of the quasi-classical Eilen- 
berger equations in type II superconductors with strong 
paramagnetic contribution. 

In this Letter we present an analytic derivation of the 
magnetic field distribution and VL form factor taking 
into account the paramagnetic effects. Our analysis is 
based on the electrodynamic theory of the Abrikosov lat- 
tice in the superconducting state where the diamagnetic 
superfluid currents are mainly determined by the Zeeman 
interaction of the electron spins with the local magnetic 
field as developed in the papeJ^^. 

The orbital and the Zeeman respective contribution is 
quantified by the Maki parameter um = '\/2i?c2o/^^p, 
where -ffc20 — '/'o/27r^o is the orbital critical field while 
ffp — Ao/a/2m is the Pauli limiting field, (/iq — 2.07 x 
10~^ G ■ err? is the fiux quantum. Unlike the majority 
of superconductors, the Pauli limiting field in CeCoIns is 
smaller'^ than the orbital critical field by a factor of > 3. 
Hence, the Zeeman interaction plays an important role 
in the mixed state field and current distributions. We 
demonstrate analytically how the VL form factor, which 
decreases with increasing magnetic field in the high tem- 



peratures region of the phase diagram, at lower temper- 
atures turns to the behaviour increasing with increasing 
field. 

We shall consider a square VL with wave vector q — 
B / (f)Q^ formed in a tetragonal type II superconduc- 
tor under magnetic field directed along the c-axis. The 
magnetic induction B = h is determined as the spatial 
average of the local magnetic field h = V x A. The GL 
theory for the VL form factor, valid in the limit k ^ 1 
and for the external field not too close to iJc2 , was devel- 
oped by J. Clem.SI Starting from the general form of the 
order parameter for an isolated vortex 



A(r) = Aoo/(r)e 



(1) 



[if is the angle measured from one of the axes in the ab 
plane), he proposed to model /(r) by the trial function 



f{r) = 



i?' 



(2) 



with R — ^7-2 -|- ^2 ^ The variational parameter was 
constrained to minimize the vortex total energy and was 
found to be in the large k limit = •\/2^, where ^ is the 
coherence length defined below. He has calculated the 
field distribution due to the orbital current and obtained 
the form factor 



B 
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where Q = \/ -|- Ki{x) is the modified Bessel 
function of first order,E21 and A is the London penetra- 
tion depth. One can write an approximative form of it 
in the conditions k ^ 1, q 3> A^^, and g^i, <C 1, 
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ln7 = C ~ 0.577 is the Euler constant. 

The found form factor, that quite slowly decreases with 
magnetic field, is reliable if one neglects the paramagnetic 
interaction of the electron spins with the magnetic field. 
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The latter leads to two extra features that are important 
in the case of a large enough Maki parameter. First, the 
two characteristic lengths ^ and A in the above expres- 
sion proves to be magnetic field dependent. Second, a 
new mechanism originating from the Zeeman interaction 
gives rise to the main contribution to the diamagnetic 
screenin^i^l the high magnetic field region of the phase 
diagram. To find it we start with the Ginzburg-Landau 
formulation including the paramagnetic effects. 

The superconductor CeCoIns has pairing symmetry 
dj,2_j,2,^^ with order parameter 

Ak(r) = ?A(k)A(r), V^(k) = ^cos(2^). (5) 

The free energy of the system is given by the Ginzburg- 
Landau (GL) functional 



J" = / d"r 



{a + e{h, - B))\Af 
-l-/3|A|4 + 7|DApy 
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where D = — iV + 2eA is the gauge-invariant derivative 
(from here we put h — c — I), and the coefficients are 
functions of temp erature and induction. In the clean 
limit they ar^i^EI! 

a = A^o(ln(r/rc) -f 5Re^'(w) - ^'(1/2)), 
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where is the digamma function, ^'^'"^(w) are its 

derivatives called by the polygamma functions, and 
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The coefficient e is proportional to B. Hence, the cor- 
responding term in the functional is negligibly small in 
the ordinary GL region near Tc (i? — > 0). In a quasi- 
two-dimensional case we deal in the first approximation 
with a cylindrical Fermi surface. For the order param- 
eter given above the average over the Fermi surface is 
(IV'(k)l^) = 3/2. 

The stationary condition for the functional in respect 
of the order parameter gives the non-linear GL equation 



aiA + 2/3|ApA + -fT>^A = 0. 



(7) 



The coherence length is defined as ^ = ^/jj\ai\, with 
ai ~ a + e{hz — B). This expression is meaningful away 
from the region where 7 = (close to which the FFLO 
phase may exist). 

The variational solution of the GL equation for an iso- 
lated vortejP is given by eqns , ^ with the order pa- 
rameter amplitude Aqo = -y/jai |/2/3. In the high n limit. 



we have h B almost everywhere. Therefore as far as 
the configuration of the order parameter is concerned, we 
will take ai ~ a. 

The stationary condition with respect to the vector 
potential gives the Maxwell equation 



1 
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V X h = jorh 
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For an isolated vortex with order parameter given by 
eqns. ([T]), ([2]) the vector potential is of the form A„(r) = 
Ay(r)(^, the orbital density of current is 



iorh 



Be27(A,(r) 



27rr • 



while the Zeeman current^^ found from ^ is 

ar 
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(10) 



Hence, we come to the equation that determines the vec- 
tor potential A„(r) 
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where A = -^Z (3 / 16we'^'y\a\ is the penetration depth. 
Writing (11) for the shifted potential a{r) = Ay{r) — 
(/)o/27rr, we obtain the differential equation with an in- 
homogeneous term of Zeeman origin 



d A d , p ^ ,2 

— (- — (ra)j - --jO = -47reA^ 
dr^rdr ' A^ ^ 
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The general solution of this equation a(r) = afi{r)+ai{r) 
consists of the sum of the solutions of the homogeneous 
and the inhomogeneous equations. The former deter- 
mines the orbital part Aort — '/'o/2t -1- of the vector 
potential Ay — Aorb + Az , 
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such that the corresponding magnetic field ho^b = ho^bZ 
iP 
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and the form factor is determined by eqn.(|3]). Taking 
into account the expressions for a, /3 and 7 we obtain 
field dependence for the orbital part of the form factor 
determined at k 1 by the first term in Q 
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The Zeeman part of the vector potential is given by 
Az{r) = a,(r) = -Ki{R/\)C{R/X) 



(16) 



with 



(21) and (15) is 
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The corresponding magnetic field = hzz takes the 
form 

hz = \ [Ko{R/X)CiR/X) + Ki{R/X)C'{R/X)] . (17) 
A 

In order to find a simpler expression for this new term, we 
can look at its behaviour in the region not too far from 
the vortex core i.e. for R X, corresponding to the 
region between the vortices of the lattice. By using the 
Bessel functions asymptotic expressions in that limit 
we find the dominating term in the eqn. (17) 



(18) 



We see that hz{r) is concentrated near the core of a 
vortex, the characteristic length associated to it be- 
ing of the order of ^y. As a remark, this term may 
also be derived by considering the equation V x = 
4:TTjz, that is valid in the absence of the orbital cur- 
rent. From ([l7| one can find the correction to (18 1, 
Shz = 47reAj^o(-R/A) ln(i?/^i,)/K^, which is small in 
the high k limit. Let us therefore consider expression 
( 18 ) for deriving the new form factor. 



The Fourier transform of the magnetic field around a 
single vortex is 



hz{q)^ J d^rhz{r)e-'^-''^87T^eAl,eyKo{qi.). 



(19) 



Hence, the contribution to the form factor that originates 
from the interaction of the electron spins with the local 
magnetic field for an array of B/(f>Q vortices per cm^ is 

Fz = 47T^^Al,^m^'{w)CKo{q^,) (20) 

One can evaluate the total energy of the magnetic field 
(quantities integrated over the whole 2-D plane) and de- 
duce that the energy of a single vortex is (in the high k 
limit) dominated by the term horb- Therefore the mini- 
mization of the energy of a single vortex gives the same 
variational parameter = as exposed in the Clem 
paper.'^ As a result we obtain 

Fz = '-^^^^m^'MKo(V2qO (21) 
For ^ 1 and < 2ttT the ratio of two form factors 
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here C,{x) is the Ricmann zeta function. We observe that 
the Zeeman part of the form factor prevails over its or- 
bital part in the phase diagram region where ^B ~ T. 

We are now able to analyze the behaviour of the total 
form factor 



F — Forb + F- 



(23) 



where Forb and Fz are given by equations (|3]) and (21) 
correspondingly. It is drawn in Fig. [T] In numerical 
calculations we assumed the values /x = gus/^ = ^J■B 
for the electron magnetic moment in the material, and 
vf — 5 X YQ^cm/s for the Fermi velocity inside the 
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FIG. 1: (Above) CeCoIns phase diagram for _ff||c-axis. The 
color lines represent the temperatures where we applied the 
model. (Below) Variations of the squared form factor at 
different temperatures including both the orbital and Zeeman 
contributions. The dashed lines represent the variations of the 
orbital part alone. 
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superconducting phase. A value for vp slightly bigger 
was given irP^ as a result of measurements of the upper 
critical field Hc2 near T^.. The 2-D density of states on 
the Fermi surface is independent of vp and is given by 
Nq = m* /2TTic, where we considered m* = 100 me for 
the electron effective mass, and £c = 7.6 x 10^^ cm is the 
lattice c-axis spacing. 

In the form factor variations, there is first a domination 
of the orbital part in the low magnetic field region {Fz 
vanishes at B = 0). We observe next a crossover to a 
region where the paramagnetic term is dominant. The 
regime where ^ goes exactly against (21| is likely to 
explain the observed constant logarithm of the squared 
form factor'^ in the interval B — 0.5 — 2T. In addition, 
these features were observed in an experiment realised 
on the s-wave superconductor TmNi2 620.1^ At larger 
fields and smaller temperatures region, the form factor 
increases with field. 

The observed form factoi^^Il falling towards zero near 
the phase transition line, where the field variation am- 
plitude decreases with the order parameter, is out of the 
region of applicability of the developed theory where iso- 
lated vortices are supposed. To describe the form factor 



behaviour reaching the maximum and decreasing at the 
approach of the upper critical field, one must match the 
present theory to the Ginzburg-Landau description valid 
in the vicinity of Hc2 at temperatures above and also 
slightly below the tricritical point where the order pa- 
rameter takes a finite value at the critical field. This 
will be published elsewhere. 

In conclusion, making use of the generalized Clem ap- 
proach we have calculated the magnetic field dependence 
of the vortex lattice form factor. Some results found 
previously in the framework of the Ginzburg-Landau- 
Abrikosov theorjffi] lead us to find a new term originat- 
ing from the interaction of the electron spins with the 
magnetic field existing inside the sample. The latter has 
a dominant contribution in the expression of the vortex 
lattice form factor at high magnetic field in the super- 
conductors with a small enough Fermi velocity.'^ The 
magnetic field caused by the diamagnetic currents origi- 
nating from the Zeeman interaction was found to be con- 
centrated in the core of vortices and gives rise to new 
features of the form factor that accounts for the mea- 
surements that are currently being made on CeCoIns.EEl 



^ A. D. Bianchi, M. Kenzelmann, L. DeBeer-Schmitt, J. 
S. White, E. M. Forgan, J. Mesot, M. Zolliker, J. 
Kohlbrecher, R. Movshovich, E. D. Bauer, J. L. Sarrao, 
Z. Fisk, C. Petrovic, M. R. Eskildsen, Science 319, 177 
(2008). 

^ J. S. White, P. Das, M. R. Eskildsen et al., arXiv: 

1001.2142 [cond-mat.supr-con]. 
^ L. DeBeer-Schmitt, C. W. Dewhurst, B. W. Hoogenboom, 

C. Petrovic, and M. R. Eskildsen, Phys. Rev. Lett. 97, 

127001 (2006). 

* J. R. Clem, J. Low. Temp. Phys. 18, 427 (1975). 

^ C. Petrovic, P. G. Pagliuso, M. F. Hundley, R. Movshovich, 
J. L. Sarrao, J. D. Thompson, Z. Fisk, and P. Monthoux, 
J. Phys.: Condens. Matter, bf 13, L337, (2001). 

® K. Izawa, H. Yamaguchi, Y. Matsuda, H. Shishido, R. Set- 
tai, and Y. Onuki, Phys. Rev. Lett. 87, 057002 2001. 
H. Aoki, T. Sakakibara, H. Shishido, R. Settai, Y. Onuki, 
P. Miranovic, and K. Machida, J. Phys.: Condens. Matter, 
bf 16, L13, (2004). 

* A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, and 
J. L. Sarrao, Phys. Rev. Lett. 91, 187004 (2003). 

^ A. Bianchi, R. Movshovich, N. Oeschler, P. Gegenwart, F. 
Steglich, J. D. Thompson, P. G. Pagliuso, and J. L. Sarrao, 
Phys. Rev. Lett. 89, 137002 (2002). 



^° C. F. Miclea, M. Nicklas, D. Parker, K. Maki, J. L. Sarrao, 
J. D. Thompson, G. Sparn, and F. Steglich, Phys. Rev. 
Lett. 96, 117001 (2006). 

" M. Ichioka and K. Machida, J. Phys.: Conf. Ser. 150, 
052074 (2009). 

M. Houzet and V. P. Mineev, Phys. Rev. B 76, 224508 
(2007). 

"Handbook of Mathematical Functions", ed. by M. 
Abramovitz and I. Stegun, Dover Publ., New York, 1965. 
^'^ A. Vorontsov and I. Vekhter, Phys. Rev. Lett. 96, 237001 
(2006). 

M. Houzet and V. P. Mineev, Phys. Rev. B 74, 144522 
(2006). 

L. DeBeer-Schmitt, M. R. Eskildsen, M. Ichioka, K. 
Machida, N. Jenkins, C. D. Dewhurst, A. B. Abraham- 
sen, S. L. Bud'ko, and P. C. Canfield, Phys. Rev. Lett. 99, 
167001 (2007). 

It is possible that at temperature close to zero, the presence 
of the quantum critical point (see J. G. Donath, F. StegUch, 
E. D. Bauer, J. L. Sarrao, and P. Gegenwart, Phys. Rev. 
Lett. 100, 136401 (2008)) enhances the effects found here 
by increasing the value of the electron effective mass. 



